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Abstract 

We review the experimental status of the 6-quark production at the Fermilab Tevatron. We 
compare all available measurements to perturbative QCD predictions (NLO and FONLL) and 
also to the parton-level cross section evaluated with parton-shower Monte Carlo generators. We 
examine both the single b cross section and the so called bb correlations. The review shows that the 
experimental situation is quite complicated because the measurements appear to be inconsistent 
among themselves. In this situation, there is no solid basis to either claim that perturbative QCD 
is challenged by these measurements or, in contrast, that long-standing discrepancies between data 
and theory have been resolved by incrementally improving the measurements and the theoretical 
prediction. 
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I. INTRODUCTION 



The bottom quark production at the Fermilab Tevatron has been called one of the few 
instances in which experimental results appear to challenge the ability of perturbative QCD 
to accurately predict absolute rates in hadronic collisions. In general, the data are underes- 
timated by the exact next-to-leading-order (NLO) QCD prediction [1, 2]. The most recent 
measurement from the Tevatron [3] is however in very good agreement with an improved 
QCD calculation (FONLL [4]), and has prompted a number of studies [5-7] suggesting that 
the apparent discrepancy has been resolved with incremental improvements of the measure- 
ments and predictions. 

Because of the experimental difficulty inherent to each result, in Sec. II we review all 
measurements of the single b cross section performed at the Tevatron, and then compare 
their average to the standard and to the improved QCD predictions. In Sec. Ill we review a 
number of measurements that compare the heavy flavor production at the Tevatron to the 
prediction of parton-shower Monte Carlo generators [8, 9]. Section IV compares cross sec- 
tions for producing both b and b quarks - centrally and above a given transverse momentum 
cut - to theoretical predictions. Our conclusions are presented in Sec. V. 

II. SINGLE B-QUARK PRODUCTION CROSS SECTION AT THE TEVATRON 

The single 6-quark cross section is inferred from the measurement of the production rate 
as a function of the transverse momentum, px, of: B hadrons; or some of their decay 
products (leptons or ip mesons); or jets produced by the hadronization of b quarks. Most 
of the Tevatron measurements correspond to b quarks produced centrally (rapidity \y b \ < 
1) and with p T > 6 GeV/c (up to p T ~ 100 GeV/c). The measured cross sections are 
tabulated as a function of the transverse momentum of different objects such as the parent 
b quark, the B hadron, or a 5-hadron decay prong (a lepton or a ip meson). This makes 
the comparison of different measurements quite difficult, and usually only a few of them 
are presented together in review articles [6, 7, 10] and compared to the same theoretical 
prediction. Therefore, we start this review with a consistency check of all available data. For 
that purpose, we use the value of the observed cross section integrated from the p T threshold 
of each experiment. We determine the ratio R of each measurement to the same theoretical 
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prediction. We then evaluate the average R and its dispersion. As standard theoretical 
prediction, we use the exact NLO calculation [1] implemented with a non-perturbative model 
for the b fragmentation 1 . The fragmentation model is based on the Peterson fragmentation 
function [11] with the e parameter set to 0.006 according to fits to e + e~ data [12]. When 
appropriate, B hadrons are decayed using the QQ Monte Carlo generator program [13]. 
The standard prediction is based on old but consistent sets of parton distribution functions 
(PDF) [MRSD [14] and MRSA [15]] since they have been used in most published works. 
There are 10 measurements of the single b cross section performed by the CDF and DjZ) 
collaborations at the Tevatron 2 . We evaluate the ratio R to the standard theory for the few 
cases in which it is not provided in the publication. The b cross sections derived from the 
production and decay of B hadrons depend on the value of f u , the fragmentation fraction, 
and the branching fractions of the B decay available at the time of publication. The value of 
these parameters has changed appreciably over time; we use the same parameters, the value 
of which will be specified in the following, for all measurements and correct accordingly the 
published cross sections. The measurements are based upon different 6-quark signatures: 

The measurement in Ref. [19] uses B mesons reconstructed through the decay B — > 
J/ijj K with J/xjj — > fi + fi~. The ratio R = 3.5 ± 15% 3 for b quarks with p™ n > 9 GeV/c is 
derived using a fragmentation fraction f u = 0.375 and a branching fraction of 5.88 x 10~ 5 [21]. 

Reference [22] is an earlier CDF measurement that uses the same decay mode and the 
same kinematic selection. Using the same fragmentation and branching fractions of the 
previous measurement, we derive R = 2.9 ± 23% 4 . 

Reference [23] presents a measurement based on the process pp — > fiX. The contribution 
of misidentified muons and of c quarks is evaluated using Monte Carlo simulations. Using 

1 We use a 6-quark mass nib — 4.75 GeV/c 2 and renormalization and factorization scales fiR = = 

2 We do not include the measurements in Refs. [16, 17] because they are based on a handful of events. The 
measurement in Ref. [18] is not included because prompt tp mesons are not separated from those produced 
by 6-quark decays. 

3 The paper quotes a discrepancy of 2.9 ± 15% with respect to the standard theory that uses the MRST [20] 
set of parton distribution functions. This discrepancy is evaluated by fitting the ratio of the data to the 
standard theory as a function of the B px- This procedure underestimates the ratio of the observed to 
predicted integrated cross sections that is 3.2 ± 15%; this ratio becomes R = 3.5 ± 15% when using the 
MRSD set of structure functions as in the measurement described next. 

4 This ratio is larger than that quoted in the publication (1.9 ± 15%) and derived by fitting the ratio of 
the data to the standard theory as a function of the B pr- 



3 



a 5-quark semileptonic branching fraction of 11.2%, the measurement yields it! = 2.5 ± 26% 
for b quarks with pf n > 21 GeV/c; for pf n > 29 GeV/c, the ratio is R = 1.9 ± 35% 5 . 

Reference [25] reports two complementary measurements that use the processes pp — > 
e D°X, with D° — > fC + 7r~, and — > e X; while the first channel is almost background 
free, the second has large background contributions of misidentified electrons and of elec- 
trons from c-quark decays. These background contributions are determined by studying the 
distribution of p^? 1 , the transverse momentum of the electron with respect to the direction 
of the momentum of all tracks around the electron direction. Using a 11.2% semileptonic 
branching fraction, the inclusive electron channel yields R = 2.4 ± 27% for b quarks with 
pmm y GeV/c. Using a branching fraction of 3.14 x 10~ 3 , the e D° channel yields 
R = 2.1 ± 34% for b quarks with pf n > 19 GeV/c 6 . 

The study in Ref. [26] uses the decay B — > J '/ijj X with J/i/j — > The B contribution 

is separated from prompt J/ip production by studying lifetime distributions. By using the 
fragmentation fraction f u = 0.375 and a branching fraction of 6.74 x 10~ 4 , the measurement 
yields R = 2.0±10% for b quarks with pf n > 9 GeV/c; for pf a > 14 GeV/c, R decreases 
to 1.7 ± 15%. 

Reference [3] reports the first measurement at y^s = 1.96 TeV 7 through the decay 
B — > J/ip X with J/ip — > This measurement extends the differential cross section 

to p T ~ GeV/c, and the data are compared only to an improved QCD calculation [5]. 
We compare to the standard theory using the information that the observed cross section is 
85% of that reported in Ref. [26], whereas it should have been 10% larger [3]. We derive a 
ratio R = 1.5 ± 9% for b quarks with pf n > 9 GeV/c; for pf n > 14 GeV/c, R decreases to 
1.3 ±9%. 

The study in Ref. [27] uses the channel pp — + n X. The b contribution is separated from 
backgrounds due to misidentified muons or c-quark decays by looking at the distribution 
of p^ 1 , the transverse momentum of the muon with respect to the direction of a jet with 
E T > 8 GeV that includes the muon. Using a semileptonic branching fraction of 11.2% and 

5 The published result uses the DFLM fits [24] . This old set of parton distribution functions is quite similar 
to the most recent PDF fits and yields theoretical cross sections that are 17% (for p™ m > 21 GeV/c) and 
11% (for p^ in > 29 GeV/c) higher than those obtained using the MRSD fits. 

6 The published ratios are based on the use the DFLM fits. We correct for the fact that these fits yield a 
theoretical prediction that is 22% and 18% larger, respectively, than that based on the MRSD fits. 

7 All other measurements considered in this review are performed at y/s — 1.8 TeV. 
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the MRSD set of parton distribution functions, the measurement yields it! = 2.1 ± 27% for 
b quarks with pf a > 6 GeV/c; the ratio is 1.7 ± 30% for pf a > 12 GeV/c. 

Reference [28] is a repetition of the previous measurement that uses slightly different 
kinematic cuts and an improved simulation of the b hadronization and decay. The publication 
uses the MRSR2 fits [29]. We correct for the fact that the theoretical cross sections are 
36% and 18% higher than those obtained using the MRSD fits for p™ n > 6 GeV/c and 
pmm > yi GeV/c, respectively. The measurement yields R — 2.5 ± 25% for b quarks with 
pf n > 6 GeV/c and R = 3.5 ± 25% for b quark with pf n > 12 GeV/c. 

Reference [30] compares the production of central b jets to the prediction of the standard 
theory. The measurement requires the presence of a muon within the jets and uses its p^ 1 
distribution to separate the 6-quark contribution from the background. The publication uses 
the MRSA' fits, and reports R = 2.4 ± 20% for b quarks with p^' m > 20 GeV/c ; the ratio 
decreases to R ~ 2.0 ± 30% for pf n > 40 GeV/c. 

The ratios of the data to the standard theory are summarized in Table I. Using the 
measurements listed in Table I, we derive an average ratio of the data to the standard 
theory that is < R >= 2.39; the RMS deviation of the 10 measurements in Table I is 0.54 
that in turn yields a 0.17 error on < R >. Before comparing the data to the improved QCD 
calculation, the summary of the experimental situation in Table I prompts a few additional 
remarks. 

The 0.54 RMS deviation is much larger than the measurement uncertainties (these un- 
certainties are dominated by systematic errors that are generally quoted as conservative 
estimates). When not using the four measurements based on detection of J ftp mesons, the 
average ratio becomes < R >= 2.33 with a 0.19 RMS deviation that, as expected, is smaller 
than the measurement uncertainties. The remaining measurements based on detection of 
J/ip mesons (first, second, sixth, and seventh line in Table I) yield < R >= 2.5 with a RMS 
deviation of 0.9. These four measurements are experimentally the cleanest and easiest to 
perform, and have the smallest systematic errors; however, they appear to be inconsistent 
among themselves. For the latter reason, it does not seem judicious to use only these four 
measurements as benchmarks of theoretical progresses [5, 31]. Additional data by CDF and 
DjZ) are certainly needed to clarify this situation. 

In most, but not all, measurements the shape of the observed transverse momentum dis- 
tribution is different from that of the standard theory (see values in parentheses in Table I). 
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TABLE I: Ratio R of measured single b cross sections to a prediction based on the exact NLO 
calculation (see text). The cross sections are for producing b quarks above a given transverse 
momentum p™ in . The ratios in parentheses highlight those cases in which data and theory appear 
to have different transverse momentum distributions and are not used in deriving < R >. Each 
measurement covers 6-quark momenta as large as 4 — 5 times the p^ m threshold. The measurement 
in the seventh row also covers small transverse momenta down to px — GeV/c. 
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(Dp [28]) 


2.5 ± 25% 




(3.5) 






b jets (//) 


(Dp [30]) 








2.4 ± 20% 


(2.0) 



In general, data and theory tend to agree better with increasing p?; in one case (ninth line 
of Table I) they disagree more. It could be a real effect, but it remains open the possibility 
that some measurements do not model correctly the background contribution as a function 
of the 6-quark transverse momentum. 

As noted in Refs. [5, 6], the measurement with b jets, listed in the last row of Table I, 
depends little on the modeling of the 6-quark fragmentation. This measurement yields 
R = 2.4 ±20%, whereas < R >= 2.39, and does not provide, in contrast to what claimed in 
Refs. [6, 7], any evidence that the b fragmentation function is a major cause of discrepancy 
between data and standard theory 8 . 

8 This result is confirmed by a recent measurement [32] of the single b cross that uses central jets with 
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The NLO prediction depends strongly on the choice of the factorization and normalization 
scales; by changing the scales by a factor of two the prediction changes by approximately 
40% [1, 2, 10]. At perturbation level, the large scale dependence of the NLO prediction 
is generally taken as a symptom of large higher-order corrections 9 . In addition, there are 
logarithmic corrections that are present at all orders of perturbation theory [33-36]. The re- 
summation of the logarithms of (pr/n^b) with next-to-leading logarithmic accuracy (NLL), 
and the matching with the fixed-order NLO calculation (FONLL) for massive quark has 
been performed in Ref. [4]. A calculation with the same level of accuracy, available for the 
production of b quarks at e + e~ colliders [37], has been used to extract non-perturbative 
fragmentation functions from LEP and SLC data [38]. These new fragmentation functions 
have been consistently 10 convoluted with the FONLL calculation to predict the B cross 
section at the Tevatron [31]. The inclusion of NLL logarithms has a modest effect in the 
Pt range considered in this review; the new fragmentation functions are harder than the 
Peterson fragmentation function and explain most of the 30% increase of the FONLL pre- 
diction with respect to the standard theory [31]. In the pt range considered in this study, 
the latest PDF fits, that include HERA data [39, 40] and a more accurate value of a s at the 
Z mass, increase by 20% the predicted 6-quark cross section with respect to the PDF fits 
used for the comparisons in Table I [6]. By also using /„ = 0.39 in place of 0.375, the final 
FONLL prediction is approximately 60% higher than the standard NLO prediction. In con- 
clusion, the ratio of the average single b cross section measured at the Tevatron with respect 
to the FONLL prediction is approximately 1.5. The uncertainty of the FONLL prediction 
is estimated to be approximately 40% [5] 11 . Therefore, the average single b cross section 
measured at the Tevatron is within the range of values predicted by the FONLL calculation. 

transverse energy Et > 40 GeV; b jets are selected identifying the presence of displaced secondary vertices. 
This study finds that the ratio of the observed cross section to that predicted by the pythia Monte Carlo 
generator is 1.2 ± 20%; this implies that the ratio of these data to the standard NLO prediction is 
approximately 2.3 (see the discussion in the next section). 
9 It is well known that there are new partonic processes that appear first at NLO, such as gluons branching 
into b and b quarks (gluon splitting) in the final or initial state of the hard scattering. 

10 As correctly noted in Ref. [31], the Peterson form of the fragmentation function used in the standard NLO 
calculation has been tuned in conjunction with a parton-level cross section evaluated with the leading-log 
(LL) approximation of parton-shower Monte Carlo programs and should not be convoluted with a NLO 
prediction. 

11 The uncertainty is estimated by changing the normalization and factorization scales by a factor of two 
(±35%) and m b by ±0.25 GeV/c 2 (±16%). 
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Exact NLO predictions do not easily allow the full simulation of events produced at the 
Tevatron. Therefore, studies that involve 6-quark production such as top quark studies or 
searches for new physics, make use of parton-shower Monte Carlo programs [8, 9]. Parton- 
level cross sections, evaluated by these generators using the leading-log (LL) approximation, 
also have large uncertainties, comparable to that of the NLO or FONLL prediction: gluon 
splitting to heavy quarks in the final state has a 30% uncertainty [41] whereas gluon splitting 
in the initial state (flavor excitation diagrams) depends on the PDF choice and can vary by 
as much as ±40% when using a wide range of structure functions in the PDF library [42]. 
Since, as correctly noted in Ref. [6], studies searching for new physics cannot depend on 
the prediction of a QCD calculation, a significant effort has been put in calibrating the 
parton-level cross section predicted by parton-shower Monte Carlo programs by using jet 
data [43-45]. Buried in top quark studies or dubious hints of new physics, the significance 
of this calibration has been overlooked. Therefore, we review it in detail in the next section. 

III. COMPARISONS WITH THE HERWIG AND PYTHIA PREDICTIONS 

It was first reported in Ref. [46] that parton-shower Monte Carlos, such as the pythia 
and HERWIG generators, predict a parton-level single b cross section that approximately 
matches the Tevatron measurements for b quarks with p? > 6 GeV/c and \y\ < 1. The 
parton-level cross section estimated with LL generators is approximately a factor of two 
larger than the exact NLO prediction because the contribution of terms of order higher than 
a 2 s is a factor of two larger than the contribution of o? s terms estimated with the exact NLO 
calculation [2]. 

Leading-order (LO) and higher-than-LO terms are sources of b and b quarks with quite 
different topological structure. The production of events with both a b and b quark with 
Pt > 6 GeV/c and \y\ < 1 is dominated by LO diagrams and the parton-level cross sections 
predicted by the exact NLO calculation is comparable to that predicted by LL Monte Carlo 
generators. At the time, both LL and NLO predictions underestimated by a factor of two the 
available measurements [28, 47, 48] 12 . Therefore, the fact that LL generators model correctly 
the single b cross section was considered merely accidental, and the source of the discrepancy 

12 These measurement identify b quarks through their semileptonic decays. 
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between data and NLO prediction was searched in non-perturbative fragmentation effects 
that enhance equally LO and NLO terms. 

In Refs. [43, 45] , the heavy flavor cross section evaluated with the herwig Monte Carlo 
generator has been tuned by using jet data collected by the CDF experiment at the Tevatron. 
This study uses four samples of data consisting of events with two or more jets, one of which 
is central and has transverse energy Et larger than 20, 50, 70, and 100 GeV, respectively, 
and a data sample, richer in heavy flavor, collected requiring two or more central jets with 
E T > 15 GeV, one of which contains a lepton with p T > 8 GeV/c from heavy-flavor decays. 
Jets containing heavy flavor are identified by finding displaced secondary vertices produced 
by the decay of b and c hadrons inside a jet; an additional algorithm uses track impact 
parameters to select jets with a small probability of originating from the primary event 
vertex. In the data, the b- and c-quark contributions are separated because both algorithms 
have the same tagging efficiency for b jets, whereas for c jets the efficiency of the second 
algorithm is approximately 2.5 times larger than that of the first algorithm. The tagging 
rates in the data are compared to those of samples simulated using the herwig Monte Carlo 
program 13 . The study compares momentum distributions of leptons or of the system of 
tracks forming a secondary vertex (decay products of the B hadron inside the jet) in the 
data and simulation. This comparison shows that the hadronization of heavy quarks at the 
Tevatron is modeled correctly by HERWIG tuned with e + e~ data. Therefore, one is allowed 
to tune the parton-level cross section predicted by the Monte Carlo generator to match the 
heavy-flavor content - or the tagging rate - of the data. The contribution of LO and higher- 
order terms can be separated because the 90% of the LO contribution consists of events 
which contain two jets with heavy flavor inside the kinematic cuts. In contrast, only 10% 
of the events due to higher-order terms contains two jets with heavy flavor in the detector 
acceptance. The higher-order contributions due to gluons splitting into heavy quarks in the 
initial and final state are disentangled by studying the 5R = ^ (S(f>) 2 + (<5?7) 2 distribution 
between two jets with heavy flavor (gluon splitting in the final state clusters at small values 
of 8R). 

References [43, 45] show that the data can be modeled by tuning the parton-level cross 
section predicted by the herwig generator within the theoretical and experimental uncer- 

13 The study uses option 1500 of version 5.6 with the MRS(G) set of parton distribution functions [15]. 
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tainties 14 . In the tuned LL generator, the contribution of higher-order terms to the single 
b cross section is approximately four times larger than the LO contribution. In contrast, for 
the same kinematics, the exact NLO calculation with standard scales returns contribu- 
tions that are only a factor of two larger than the a 2 s contribution. As discussed in the next 
section, the study of bb correlations can be used to assess the correct ratio of higher-than-LO 
to LO contributions. 

IV. MEASUREMENT OF THE BB CORRELATIONS 

The cross section for producing both b and b quarks centrally and above a given pt 
threshold, a bb or bb correlation, is dominated by LO terms, and the LL and NLO predictions 
are quite close 15 . In addition, the exact NLO prediction of a bb depends little on the choice 
of the normalization and factorization scales as well as on the 6-quark mass 16 and appears 
to be a robust prediction of perturbative QCD. 

Therefore, it is important to determine precisely the value of R 2 b, the ratio of a bb measured 
at the Tevatron to the exact NLO prediction (or to the LL prediction that is very close). 
A ratio R 2 b — 1 would imply that the parton-level cross section predicted by LL Monte 
Carlo generators is correct and that the contribution of higher-than-LO terms has to be a 
factor of two larger than in the present NLO or FONLL prediction. If the ratio R 2b is much 
larger than one, then the agreement between the observed single b cross section and the 
prediction of LL Monte Carlo generators is fortuitous. Since the NLO prediction of a bb is 
robust, agreement with the data may be found by using harder fragmentation functions as in 
the FONLL calculation. Unfortunately, the status of the a bb measurements at the Tevatron 
is quite disconcerting. 

The study in Ref. [43] (CDF) uses two central jets with Et > 15 GeV, each containing a 
secondary vertex due to b- or 6-quark decays. The LL prediction, tuned to fit the data, yields 

14 The gluon splitting in the final state predicted by herwig has to be increased by (40 ± 20)%. Before 
tuning the simulation, the size of gluon splitting in the final state predicted by herwig is 1/2 of that in 
the initial state. 

15 For example, in Ref [43] the tuned LL generator predicts a modest contribution of higher-than-LO order 
terms to <j bb (~ 30%); for the same kinematics, the exact NLO calculation predicts a ~ 15% contribution 
of higher-than-LO order terms. In this case, the LL and NLO predictions of a bb are within 20%. 

16 The prediction changes by no more than 15% by changing the scales by a factor of two and m b by ± 0.25 
GeV/c 2 [43]. 
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R 2b = 1.2 15 with a 25% uncertainty mostly due to the efficiency for finding a secondary 
vertex in a jet. 

A recent measurement [49] (CDF) supports the conclusion of Ref. [43]. The study in 
Ref. [49] uses events containing two central jets with > 30 and 20 GeV, respectively; 
pairs of b jets are also identified by requiring the presence of displaced secondary vertices. 
This study finds that the ratio of a bb to the LL pythia prediction is 0.9 ± 31%, while the 
ratio of a bb to the NLO prediction 17 is R 2b = 1.0 ± 32%. 

In contrast, discrepancies between data and the NLO prediction of a bb are observed when 
identifying b quarks through their semileptonic decay into muons. 

The study in Ref. [47] (CDF) uses events with a muon recoiling against a jet that contains 
tracks with large impact parameter (b jet). Using the average branching fraction BR = 10.3 
for b — > ji X decays and 10.2 for b — > c X — > /i Y sequential decays [51], the ratio of a bb to 
the exact NLO prediction is measured to be R 2b = 1-5 ± 10% for b and b quarks produced 
centrally and with p^ in > 12 GeV/c. 

Reference [48] (CDF) reports a measurement that uses two centrally produced muons. 
By using the square of the semileptonic branching fraction quoted above, the study yields 
R 2b = 3.0 ± 20% for central b and b quarks with p™ in > 6 GeV/c. 

Reference [28] (DO) reports an analogous measurement that also uses two centrally pro- 
duced muons. The square of the semileptonic branching fraction is evaluated with the 
isajet generator [52] implemented with the QQ decay table and is consistent with the value 
used by CDF. The study yields a ratio R 2b = 2.3 ± 33% for central b and b quarks with 
pf n > 7 GeV/c. 

These five measurements, listed in Table II, yield < R 2b >= 1.8 with a 0.8 RMS devia- 
tion. Such a large RMS deviation indicates that the experimental results are inconsistent 
among themselves. Additional measurements are certainly needed to clarify the experimen- 
tal situation. However, it is quite obvious that the present discrepancies are reduced if the 
rate of observed semileptonic decays is approximately 50% higher than what is expected 
because: (a) lepton identification efficiencies are underestimated by approximately 50% or 
(b) additional objects with a 100% semileptonic branching ratio and a cross section of the 
order of 1/10 of the b cross section are produced [43]. Reference [43] has investigated these 

17 In this case the NLO prediction has been evaluated using the MC@NLO Monte Carlo generator [50]. 
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TABLE II: Ratio i?26 of a b i, the observed cross section for producing both b and b quarks, centrally 
and above a given p™ in threshold, to the exact NLO prediction (see text). Each measurement covers 
6-quark momenta as large as 4 — 5 times the p™ m threshold. Jets produced by b and b quarks are 
identified by the presence of displaced secondary vertex or tracks with a large impact parameter. 
Muons from b and b decays are separated from the background by studying impact parameter [48] 
or p™ 1 [28] distributions. 



channel (experiment) i?2b for p™ m (GeV/c) > 

6-7 10 15 -20 

6 + 6 jets (CDF [43]) 1.2 ±25% 

6 + 6 jets (CDF [49]) 1.0 ± 32% 

H + b jet (CDF [47]) 1.5 ± 10% 

+ (CDF [48]) 3.0 ±20% 

^+±^- (Dp [28]) 2.3 ± 



hypotheses by comparing the rate of observed and predicted leptons from 6-quark decays in 
jets that recoil against a generic jet or a jet that also contains a lepton (the jets are central 
with E T > 15 GeV). This study finds that in the second case the rate of jets containing a 
lepton from presumed b decays is 50% higher than in the first case. The magnitude of the 
effect is consistent with hypothesis (b). In light of this observation, it is worth to go back to 
Table I. If there was a reason to disregard the first two measurements using the B — > J / ipK 
channel, one would see a completely different picture. Six measurements identify b quarks 
through their semileptonic decay and yield < R >= 2.33 ± 0.06; the measurements with 
inclusive J/tp mesons (sixth and seventh line) do not use b semileptonic decays and yield 
< R >= 1.7 ± 0.1; this conjecture very well highlights the need for additional cross-checks 
of the measurements based onB^ J/ipK decays and on the inclusive J/ip production. 



V. CONCLUSIONS 



We review all measurements of the single b cross section performed at the Tevatron and 
compare them to an exact NLO perturbative QCD prediction, that uses pre-HERA sets 
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of parton distribution functions and the Peterson fragmentation function, in order to test 
their consistency. We also compare the data to an improved QCD calculation (FONLL) 
and to the prediction of LL Monte Carlo generators. The average ratio of the data to the 
NLO prediction is < R >= 2.39 with a 0.54 RMS deviation. The RMS deviation is much 
larger than the quoted measurement uncertainties, and indicates that experimental results 
are inconsistent among themselves. With this caveat, the average of the data is found to 
be in agreement with the parton-level cross section evaluated with parton-shower Monte 
Carlo generators and is within the range of uncertainty of the FONLL prediction that in 
turn is 60% higher than the NLO prediction. The increase of the FONLL prediction with 
respect to the NLO calculation is mostly due to PDF improvements and the usage of a 
harder fragmentation function, whereas the parton-level cross section are the same in both 
predictions. On the contrary, the contribution of higher-than-LO terms returned by LL 
Monte Carlo generators fitted to the data is approximately a factor of two larger than that 
in the FONLL or NLO calculations. The measurement of a bb) the cross section for producing 
both b and b quarks centrally and above the same pr threshold, has a decisive role in assessing 
the correct parton-level cross section. In fact, the higher-than-LO contribution to a bb is quite 
modest in all theoretical approaches. Because of the use of harder fragmentation functions, 
the FONLL calculation yields a prediction appreciably larger than that of LL generators 
or NLO generators convoluted with the Peterson fragmentation function. Unfortunately, 
the experimental situation is quite disconcerting and only raises additional questions. The 
average ratio of the a bb measurements to the NLO prediction is < R 2 b >= 1.8 with a 0.8 RMS 
deviation, and suggests that these measurements are also inconsistent among themselves. 
The < i?26 > value supports the FONLL approach. However, the level of agreement between 
data and theory appears to be a function of the number of semileptonic decays used to 
identify b quarks. In this situation, it cannot be excluded that the b parton-level cross 
section is correctly described by LL Monte Carlo generators, and that measurements using 
b semileptonic decays are affected by new physics. 
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